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Abstract 

An analytical form of the quantum magnetization oscillations (de Haas-van Alphen 
effect) is derived for two- and quasi two-dimensional metals in normal and supercon- 
ducting mixed states. The theory is developed under condition fi/uj c 1 (/i is the 
chemical potential and uj c the cyclotron frequency) , which is proved to be valid for us- 
ing grand canonical ensemble in the systems of low dimensionality. Effects of impurity, 
temperature, spin-splitting and vortex lattice - in the case of superconductors of type 
II -, are taken into account. Contrary to the three dimensional case, the oscillations 
in sufficiently pure systems of low dimensionality and at sufficiently low temperatures 
are characterized by a saw-tooth wave form, which smoothened with temperature and 
concentration of impurities growth. In the normal quasi two-dimensional systems, the 
expression for the magnetization oscillations includes an extra factor expressed through 
the transfer integral between the layers. The additional damping effect due to the vortex 
lattice is found. The criterion of proximity to the upper critical field for the observation 
of de Haas-van Alphen effect in the superconducting mixed state is established. 

1 Introduction 

De Haas- van Alphen (dHvA) effect is well known to be a powerful tool to get information 
on the shape of the Fermi surface of normal metals. In contrast to three-dimensional (3D) 
metals where the experimental observations have well established base given by the Lifshitz- 
Kosevich (1956) theory, there is no unanimously acknowledged and used theory applicable to 
2D or quasi 2D layered materials. Schoenberg (1984) introduced an expression for magnetic 
oscillations in 2D metals where he put phenomenologically the effect of temperature. Effect 
of impurities is not considered, although it can be important. There are two essential features 
of dHvA effect in 2D case commonly believed differing it from its 3D counterpart : 

(i) The sharp saw-tooth like shape of dHvA signal at sufficiently low temperatures. 

(ii) Strong dependence of the chemical potential on the magnetic field, that is why the 
proper derivation has to be done in the canonical ensemble with fixed number of particles. 

To make clear this subject we note that the property (i) takes place in 2D case both 
at constant number of particles or at constant chemical potential so long we are in the low 
temperature T <C uj c and high purity region T < uj c . Here, lo c is the distance between 
Landau levels, T is the width of level. As for property (ii) its validity is directly related to 
the concentration of the charge carriers we deal with. For low enough concentrations when 
0J C ~ £f — =Jjr (n is the electron density and m* is the effective mass), the oscillations of 
the chemical potential /j, with magnetic field are important and one has to proceed with 
the calculations in canonical ensemble (see the papers (Jauregui et at 1990, Harrisonei 
al. 1996, Grigoriev and Vagner 1999, Nakano 1998)). On the contrary under the condition 
H w sp = ™ 3> lo c the oscillations of the chemical potential are small and one may work 
making use the grand canonical ensemble as in 3D case. 
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The analytical theory of dHvA effect in 2D and layered quasi 2D metals taking into 
account the effects of temperature, impurities and spin-splitting is developed in the first 
part of this paper, ft is proved that the grand canonical ensemble is well accurate for 
systems of low dimensionality in the limit /u/w c » 1. In pure case and at zero temperature 
the oscillations have the saw-tooth shape. The finite temperature as well as the finite amount 
of impurities (even at zero temperature) lead to the natural smoothening of the oscillations. 
The results of calculations (but not the calculations themselves) of magnetization oscillations 
in frame of grand canonical ensemble for quasi 2D spinless and perfect crystal model can be 
found in the paper (Nakano 1998). If do not take in mind the opposite sign of magnetization 
(it is probably due to a misprint) they coincide with ours taken in the absence of impurities 
and spin paramagnetism. The recent experimental studies (Wosnitsa et al. 2000) in several 
organic compounds demonstrate magnetization oscillations behavior similar to that should 
be in the systems with fixed chemical potential. 

In a parallel direction to the study of dHvA effect in the normal low dimensional systems, 
the possibility of doing quantum oscillations measurements in the superconducting mixed 
state has been recognized after its first manifestation by Graebner and Robbins (Graebner 
and Robbins 1976). During the last ten years, dHvA effect has already been observed both 
in mixed state and normal state of many types of superconductors like NbSe2 (Haworth et 
al. 1999), Nb 3 Sn (Harrison et al. 1994), V 3 Si (Corcoran et al. 1994), YNi 2 B 2 C (Goll et 
al. 1996, Terashima et al. 1997), CeRu 2 (Hcdo et al. 1998), UPd 2 Al 3 (Haga et al. 1999), 
URu 2 Si 2 (Ohkuni et al. 1999), k — (BEDT — TTF) 2 Cu(NCS) 2 (Sasaki et al 1998). Apart 
the low temperature 

rp „ eH c2 T c 2 

mc sf 

the condition of observation of dHvA effect in the superconducting mixed state consists of 
ultra purity of the samples 

Here I is the mean free path of quasiparticles, £ is the coherence length. Both of these 
demands mean that all enumerated above samples of materials are ultra clean and at the 
same time very strong type-II superconductors with high upper critical field value. The 
latter is provided by the higher critical temperature and the lower Fermi energy than in 
ordinary type-II superconductors. 

The main observation consists of that the frequency of oscillations in the mixed state 
remains the same but the amplitude decreases with decreasing field more rapidly than in 
the normal state. The effect becomes unobscrvable at fields several times smaller than H c2 
but still much larger than the thermodynamic critical field H c = H c2 /k. Here k is the 
Ginzburg-Landau parameter, which is of the order of 20-30 typically for observations of 
dHvA in mixed state materials. That means the distance between the vortices is of the 
order of its core diameter in the whole region of the observation. Hence the latter can be 
considered as vicinity of the upper critical field where the Abrikosov solution for the vortex 
lattice serves as a good approximation. 

Another important observation (Goll al. 1996; T. Terashima et al. 1997) is that dHvA 
effect in the superconducting mixed state persists below the upper critical field not for all 
but for the electron orbits with relatively smaller radius (or cross section area). 

Numerous theoretical studies were done to understand the effect of vortex lattice. The 
essential point here is a nondiagonality of the order parameter matrix in the Landau repre- 
sentation preventing a simple derivation of the quasiparticle energy spectrum in the mixed 
state. The intention to avoid of this problem leads to the idea to work with oversimplifycd 
BCS type spectrum (Miller and Gyorffy 1995; Dukan and Tesanovich 1995); Gvozdikov and 
Gvozdikova 1998), which is valid either in the quasiclassical region far below of H c2 where 
the dHvA effect does not take place or in ultra quantum limit with only few filled Landau 
levels. In the latter case there are however doubts as in the applicability of the BCS theory 
as in the existance of the superconductivity. 
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It should be mentioned also an attempt to develope a theory for anisotropic supercon- 
ductors with nodes in the gap on the Fermi surface far below H C 2 where a vortex core radius 
is much smaller than the distance between vortices , and one can assume the space constancy 
of the order parameter modulus (Gorkov and Schrieffer 1998). As we pointed above it does 
not corresponds to the observations region. 

The proper derivation has been proposed by K.Maki (1991) and M.Stephen (1992). There 
were established an additional attenuation of the dHvA amplitude in the superconducting 
mixed state. The order parameter has been determined self-consistently by M.Stephen 
that has allowed to point out the field interval near H C 2 where the effect is accessible for 
observation. The mentioned results were obtained in neglect by the nondiagonal elements of 
the self-energy matrix , or sort of random vortex lattice potential approximation has been 
used. 

Somethat later the numerical calculations for two-dimensional electron system has been 
performed (Norman and MacDonald 1996; Bruun, Nicos Nicopoulos and Johnson 1997). 

Another type of analitical results ( with less rapid attenuation of dHvA signal with 
decreasing of magnetic field in the superconducting region) were represented in the paper 
(Zhuravlev et al, 1997), where so called strict phase coherence approach has been applied. 
Although the latter has definitive theoretical interest it seems to be less realistic than the 
random vortex lattice potential approximation has been used in the Stephen's paper. 

The treatments (Stephen 1992, Norman and McDonald 1996, Bruun, Nicos Nicopoulos 
and Johnson 1997, Zhuravlev et al, 1997) were limited by the condition T > lu c /2it 2 which 
is necessary for the convergency at low temperatures. The spreading of the theory on the 
low temperature region to c > T > T by an introduction of the width of the Landau levels 
r originating of impurity scattering has been done in the papers (Vavilov and Mineev 1997; 
1998; Mineev 1999). It was shown that a superconducting state is gapless in the mixed state 
region below the upper critical field if 



where £f/^c is the number of Landau levels below the Fermi surface at H = H C 2- This 
field interval is negligible for any typical type-II superconductor. However for those partic- 
ular ultraclean materials with very high H C 2 and very small £f, where dHvA effect in the 
superconducting mixed state have been observed, the value of ^Jef /^>c2 is of the order of or 
even smaller than ten and the presented theory has sizable region of applicability below the 
upper critical field. In this region, the oscillatory contribution to the density of states and 
the damping of the amplitude of the magnetization oscillation in the superconducting state 
were found. 

It also important to note that amount of Landau levels below the Fermi energy is different 
for the different bands (e f in (0) should be diminished on the value of energy at the center 
of the band under consideration). That is why the condition ([!]) is less restrictive for the 
electron orbits with smaller radius or cross section area. This fact is in correspondence with 
the observations mentioned above. 

In the unconventional superconducting states where the lines of extremal cross-section 
of Fermi-surface by the plane perpendicular to magnetic field coincide with lines of zeros of 
the order parameter the amplitude of the magnetization oscillations is practically the same 
as in a normal state. In the superconducting state with an other distribution of zeros the 
damping of dHvA oscillations corresponds qualitatively to the ordinary superconductivity. 

The further development of the theory of the dHvA effect in the superconducting mixed 
state , has been done in the paper (Mineev 2000). In three-dimentional isotropic model there 
was analitically derived the Landau expansion for the free energy of the superconducting 
mixed state near the upper critical field in powers of the square modulus of the order 
parameter averaged over Abrikosov lattice with Landau quantization of the quasiparticle 
energy taken into consideration. 




(1) 
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The condition of the valitity of such the expantion were established and it is given by 
the formula (|l|). 

As in the normal metal the oscillating with magnetic field terms represent the tiny 
corrections to the nonoscillating free energy of the superconducting state. However they 
are fast oscillating functions of the ratio 2it£f/lo c . Becouse of that after the differentiation 
over the magnetic field M = — T the oscillating part of the magnetization in the 

superconducting mixed state proves to be even larger than corresponding nonoscillating 
part of magnetization. 

The analytical theory of the dHvA effect in the superconducting mixed state in the 
vicinity of the upper critical field for 2D and quasi 2D layered superconductors is developed 
in the second part of the present article. Our treatment follows the paper (Mineev 2000) 
where the corresponding calculations have been performed for 3D superconductors. 



2 De Haas- van Alphen effect in two-dimensional normal 
metals 

In order to derive the magnetization M, it is necessary to calculate first the free energy F 
which is related to the thermodynamic potential VI in grand canonical ensemble byQ 

F = fl + fiN (2) 

with 

/>oo 

n = -T 3(e) ln(l + e^)de (3) 
J o 

where y, is the chemical potential of the system, N is the number of electrons per unit volume 
and g(e) is the density of states in presence of the magnetic field H and the impurities. 



2.1 Density of States 

The latter can be achieved from the imaginary part of the space averaged single-electron 
Green's function G(e) using the relation 

g{e) = —sign(e- /x)SG(e). (4) 

We then follow the proceeding of Bychkov (1961), who has dealt with the effect of impurities 
in 3D metals under magnetic field using a perturbation method. The broadening of each 
Landau level into the energy distribution with Lorenz density of states has been found in 
this paper supposed to be not so bad aproximation in 2D and quasi 2D cases apparently 
valid in not very strong magnetic fields (see Prange 1987). So, we shall use the Green's 
function in the following form 



G ® ^ e - E m + iTsignie - n) (5) 

m " 

where m specifies the quantum numbers, L the width of Landau's levels due to impurity 
scattering (F = l/2r with r a mean free time of quasiparticles) and E m the energy spectrum 
of the stationary states in absence of impurities. The quantum numbers depend on the 
dimensionality of the considered system : for the 2D case, the set of quantum numbers m 
includes the magnetic quantum number n of the Landau's levels, a wave number k y which 
is connected with the center of the orbit of the electrons, and the spin a = ±1 ; for the 3D 

1 Planck's constant h, the Boltzmann's constant kg and the ratio of electron charge |e| to the velocity of 
light c are chosen to be unity throughout the paper 
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case, we have in addition the momentum k z . The field H is directed along the z axis, that 
is in 2D case perpendicular to the conducting plane. The energies are : 



El d kz a =uj c (n+ l -) + ^ + a^H, (6) 

E% tr = u e (n+^)+<Jii e H (7) 

One needs to remember that all the energies should be counted from the energy eb at the 
center of the considered band. Keeping in mind this circumstance we shall shift the chemical 
potential in final formulas on this value. The expression for the single Green's function for 
a given spin a becomes in the 2D case 



& (e) = N u; c J2 r—T-, 5— ^ ■ jf \ • (8) 

^ Q e — uj c (n + 2 j — cr^ e H + iL sign(e — fi) 

Here Nq = -p- is the normal metal electron density of states per one spin projection and \x e 
is the electron's magnetic moment. The expression (^) is divergent. In fact, it is the case 
only for the real part of G a {s) ; the imaginary part, which we have to calculate next, is 
convergent 

o/=w/v> at -Tsignje - fx) 

. G ( £ )=iVo. c g (£ _ cJc(n+|) _^ )2+r2 . 0) 

The Poisson's summation formula 

oo +oo r +oo 



£/w = E / /W e 

n=0 Z=-oo a 



where a is a number between -1 and 0, can now be used properly. Writing e = £ + fj,, 
changing variable t to x = t — , and then setting the lower limit of the integral equal to 
— oo since we consider the limit 3> 1, we get 

9?G CT ( £ ) = -TrN signC + ^G^ sc (e) (10) 

with 

ZG" osc (e) = -N u c Tsignt; ]T e" 2 ^ / -— ^ TTTJ— 7^ (U) 



The first term in the right hand side of equation ( jio| ) comes from the contribution of Z = in 
the sum and corresponds to the result in absence of field. The other contributions I give 
an oscillatory part depending on the field. The remaining integration is quite straightforward 
by a residue calculation and yields 

= 2nsign(e - v)N J2(-l) l+1 cos ~ a ^ H) \ e^ 1 ^ . (12) 

i=i ^ UJ ° ' 

From relation (Q), we get directly the searched density of states, after summing up the two 
spin states 



g 2d (e) = 2N 



+oo „ 

1 + 2 VY-1) 1 cos(2tt/ — ) cos(27rZ^) e" 27 ^ 

1=1 C C 



(13) 
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2.2 Chemical potential 

Making use of equation (||) , the oscillatory part of the thermodynamic potential f2 is per- 
formed by an integration by parts and by the change of variable x = ^ £ Z/^ 



St-^gfc^ .-icosp^i r (14) 

^ I Uc ^ J-J± l+e x 

The lower limit of the integral may be set equal to — oo since in our investigation we have 
the condition ji^>T. Using the value of the integral 



ay 



1 + e y sinh(a7r) 



we find at last 



^ (-1)' o i r u cos(2vr/^) 

= 2iV T Wc ^ i-1 cos 2?^ , ; "° ; . (15) 

The nonoscillatory part of the thermodynamic potential gives on its side 

fig* = -N Qf i 2 . (16) 

The number of electrons is defined by N — — (j&^J an( l thus 

+2? cos(27r/^) u 
N = 2N^ + 47riV rV(-l) i , ) %' sm(2nl-^) e" M - . (17) 

This formula determines of the number of particles as the function of (H, T) at fixed /i. On 
the other hand one can consider it as the equation for the chemical potential as the function 
of (H, T) at fixed number of particles. Recognizing the Fermi energy 

N 

£f = (18) 

one can rewright this equation as 



l=i 

with 



e F + ^Y A- sin(2^) cos( 27 r/^) e" 2 ^ (19) 

7r f— ' I sinhA; uj c lj c 



2ir 2 lT 

A; = . 



We see that at ef w c the second oscillating term represent the small correction to the 
const value of chemical potential. On the other hand the expression (|ll|) has been obtained 
in neglect of terms of the order of w c //i. That is why it would be out of limits our accuracy 
to keep the difference between the chemical potential and Fermi energy in the right hand 
side of (p~9|) . So, finally we get 

fl = £p + i ±y til * sin (2^) cos(27r^) e" 2 ^ . (20) 

7r I sinn A; uj c uj c 
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It is worth to remind here that in real crystal one should to shift all the energies on the 
energy in center of band: fi — ► [i — £(,, sf —> £f — £&• 

At zero temperature, ignoring for brevity the spin splitting and summing up over I, we 
find a simple expression 



U! 



sin(27r^) 



/i = £f + — arctan ( — — ] (21) 



7T 



v cos(2tt^) +e 



which is obviously has at Y — saw-tooth shape. At any finite Y, the oscillations of the 
chemical potential are smooth. For instance for 2ttY/ui c > 1 

V = £ F + — sin(2^— )e- 2n ^. (22) 

7T L0 C 

On this stage one should stress ones more that all the results has been obtained here at 
the condition fj,/ui c ^> 1- 



2.3 Magnetization 

The magnetization M = — (^y) T can be found from the relation (14) . Keeping only the 
preponderant terms and ignoring as before the small difference between [i and ep in the 
arguments of oscillating harmonics R, we obtain 



M% sc = $>J?' = N ^e F j: { -^-^-sH2^)cos(2nl^)e-^. (23) 
L — ' 7r H * — ' I Sinn A; oj r ui r 

i=i i=i 

This formula without impurities factor corresponds to the intuitive formula given by 
Schoenberg (1984). After ( p0| ) and ( p3| ) we notice that at ep lo c the magnetization is 
directly related to the oscillating part of /i by 

<* ac = ^/W = M ^. (24) 

ti UJ C 

This result, which is valid at any temperature smaller than the ep, generalizes that was 
given in the recent paper (Itskovsky, Maniv and Vagner 2000) for the range of temperature 

T < LO c . 

For T=0, this expression is simplified since a summation over the integers I is possible. 
Without spin-splitting, we have 



r2d 2N w c _J sin(27r^) 



MZsc = —-j?e F arctan F \ 2ff r (25) 

7T H \_ cos (27r^) + e 27T ^ J 

and with spin-splitting 

^ N w c ( sin(4 7 r^) + 2sin(27rg)co S (2 7 r^f)e 2 ^ \ 

J"r. o«,. = —£f arctan - - — = =- . 26 

' « H ^08(4^^) + 2cos(27rff) cos(27r^)e 2 "- +e 4 "-/ 



For 2ttY/lli c > 1, this yields 



MlU = ^^e,sin(2^)cos(2 7 r^ ) e- 2 ^. (27) 

7T tl UJ C U) c 



2 The restoration of the oscillating difference fi — ep would be resulted in the additional frequences in 
Fourier spectrum of the magnetization which certainly present when the ratio to/ep is not small. One should 
stress however that both formulas for the chemical potential and for the magnetization are derived here by 
making use the expression for the density of states have been found in the neglect of terms of the order of 
ld c /ef and keeping of such the terms in the oscillating arguments would be out the limits of accuracy we 
used. 
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3 De Haas-van Alphen effect in quasi two-dimensional 
normal metals 



In layered metals, the simplest spectrum for electrons is 

E n ,kz,a = u c {n + ~) - 2tcos(k z d) + afi e H. (28) 

The parameter t is the transfer integral along the z axis, d is the distance between the layers. 
Here the field H is perpendicular to the layers. More generally, in a tilted magnetic field, 
the spectrum energy is quantized as (Yamaji 1989) 

E n ,k z .a = u c cos 9 (n + — ) — 2t, J Q(kpd tan 9) cos(fc 2 d) + <jfj, e H (29) 

where Jo is the Bessel's function of zero order and kp is the momentum at the Fermi's surface. 
For convenience, we will work with spectrum (p8|). The corresponding generalization for the 
case of finite tilting is easy with simple substitution of the angle dependent parameters. 
According to (|h, the space averaged Green's function becomes henceforth 



N UJ, 



2nd ^ j - 

n=0 J 



dy: 



u c {n + 5) — cr^eH + 2£ cos^ + iTsign{e — /x) 



(30) 



Here we consider that t < or ~ lu c so that roughly speaking only the modulation of one 
Landau level is taken into account. For t 3> u> c , many Landau levels play a role and we can 
do the following approximation 

1 2t , , 1 2i., (k z d) 2 . 

In this limit, we are again in the 3D case. Next, we proceed exactly like previously using 
the Poisson's summation formula. The density of states includes then an additional Bessel 
function factor 



9 q2d (e) 



2N 



1 + 2 V (-1)' cos(27rZ — ) Jo(2ttZ— ) cos(2ttZ^) e" 2 ^ 



(31) 



with 



1 r +7r 



( 27T22 COS tp 



dtp 



and the oscillatory part of the thermodynamic potential becomes 



n& = ^^f^ t! 1 M™*) 003(2,^) e-^. (32) 



^—r l sinh(Ai) 



Thus 



Mq 2d N 2u; c +^ (-1)^ A, ,fcff, _ M i 

M l,osc=-, — TF^V ; — -cos(2tt/ e 

d 7r H '-^ I sinn A; uo c 

1=1 1 c 

x (sin(27rZ-^)J (27W — ) + - cos^-^-U^tt/— ) j (33) 
where Ji(z) = —J' Q {z) is the Bessel's function of first order. 
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Assuming that the transfer energy t is much smaller than the Fermi energy ep , we obtain 



M nosc = —, 77 M > ; . , , sm(27rt — )Jo(27r/ — )cos(27r/ )e "<= . (34) 

n ' osc d n I sum A; w c w c uj c ' v ; 

Here we remind that in reality the chemical potential fi has to be reduced on the energy 
of the band center /1 — > fi — ej ps ef — £b- Without spin-splitting and impurities factor, 
this result is the same as given in (Nakano 1998) except the opposite sign probably due to 
a misprint. The extra factor Jo(2nl^) in this expression induces a rich behavior for the 
manifestly field dependent amplitude of the harmonics. In formal limit 2t 3> u> c , 




and thus we return back to expression closed to Lifshitz-Kosevich result (Lifshits and Ko- 
sevich 1956) for the amplitude of the harmonics. For 2t <C w c , Jo(z) ~ 1 and we find again 
the 2D expression (|23|). 

4 De Haas-van Alphen effect in mixed state of two- and 
quasi two-dimensional superconductors 

In this section, we want to determine the effect of the low-dimensionality on the additional 
damping factor due to the vortex lattice. In this purpose, we follow the work (Mineev 2000) 
(and also take the same notations) where the theory of dHvA effect for 3D metals in the 
superconducting mixed state has been developed, and adapt it to 2D and quasi 2D systems. 
It should be mentioned that under 2D systems we shall imply 3D layered crystals with 
negligibly small interaction between layers, such that the mixed state in the magnetic field 
perpendicular to the layers represents Abrikosov but not Pearl vortex lattice. The treatment 
of the problem is carried out in the framework of the Gorkov's formalism for a conventional 
superconductor. 



4.1 Free energy 

The free energy density in mixed state near the upper critical field H C 2 is developed in 
powers of the square modulus of the order parameter averaged over the Abrikosov's lattice 
(Mineev 2000) 

F s - F n = aA 2 + ^A 4 (35) 

where 



and 



a = - - / e — -K 2 CR) dR (36) 
9 



P=\ J r(r 1 )/(r 2 )r(r 3 )/(r 4 )^ 4 (r 1 ,r 2 ,r3,r 4 )dr 1 dr 2 dr 3 dr 4 . (37) 



Here g is the constant of the pairing interaction, p 2 = R 2 — Z 2 , 



„ . . 1 1 a v-^ I .2ttv fx -kv ,\ 2 
/r =2 1/4 ]T exp (i y-(- + —X 

y— integer v 
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is the Abrikosov square lattice solution for the order parameter where the elementary cell 
edge a is such that a 2 — n\ 2 and A = H^ 1 / 2 is the magnetic length and A is the order 
parameter amplitude. 

The functions K 2 and K 4 have the following expressions : 

K 2 (R) = \Y, r^G CT (R,^)G- CT (R,-^), (38) 

a=±l v 

G CT ( ri ,r 2 ,w„) = cxp (i ^ ' A(s) ds^j G° ( ri - t 2 ,w v ), (39) 
^4(ri,r 2 ,r 3 ,r 4 ) = ^ £ tJ^ G° (ri, r 2 , cJ,)G- ff (r 2 , r 3 , -cJ,)G ff (r s , r 4 , w„) 

<T=±1 Z/ 

xG- ff (r 1 ,r 4 ,-^), (40) 



a),, + Tsignujv, u) v = i:T{2v + 1). (41) 

To obtain this formula, it is assumed that the magnetic field is uniform and coincides with 
external field, which is a good approximation in the vicinity of H c2 for superconductors with 
large Ginzburg-Landau parameter. 

The coefficients a and /3, which depend on the dimensionality of the system by way 
of the Green's functions, are calculated respectively in appendix A and in appendix B in 
the semi-classical approximation 1. They consist of a smooth function and a fast 

oscillating function of the magnetic field : 

a(H,T)=a(H,T) + a osc (H,T), (42) 



(3(H,T)=P(H,T)+(3 osc (H,T). (43) 
For a two-dimensional system, we find 

WW-ZZ^Q. (44, 
where H c2 (T) is the upper critical field at low temperatures averaged over the oscillations, 



a 2 d sc (H, T) = -47r 3 / 2 A -L= Y7-1)' cos(2^) cos(27r(Z + 2)^) 

• /^W c ^—^ U! c U! c 

-2tt(1+2) 

sinh A( +2 



x (45) 



and 

R 2d (H T — s ! — AT„ ,,, 

16^ (r 2 + (neH) 2 ) 2 IQkuT 2 



(3 2d (H,T = 0) = N " c i?"/ Me 5?l - ^o,^ (46) 



which is valid for fi e H < T < uo c , 

+oo p 

P 2 J C (H, T = 0) = A -^ ]T(-1)< cos( 27 r/ ^)/(2 7 r/-) e - 2 -^ (47) 

where /(a;) = J °° ^qrjp + (T+^p e_xa ^J/ i s a bounded function. Minimization of the free 
energy over A yields 

Fs=F n -f p . (48) 
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4.2 Magnetization 

Magnetization is given by M s — — Q^f- . In the differentiation we keep only the fast oscillating 
terms 

s ' osc ~ n ' osc + 8H 2(3 2 dH ■ 1 9j 

Here M n ^ osc is the oscillating part of the normal metal magnetization calculated previously 
(see equation (p3|)). Like in the 3D case (Mineev 2000), at 47iT ~ ui c the following inequality 
takes place 

a da osc a 2 d/3osc . , 

|M„, osc |>|^ w |>|-^ w | (50) 

until the new more restrictive than in the 3D case (see (|l|)) condition : 



H c2 (T) - H 



1/2 



< - ■ (51) 



H C 2o \ A 4 

With the previous condition we keep also only the first two terms and obtain for the 
oscillating part of the magnetization 

_ Q +00 

M!f osc = Ml% sc + f^f=£ M^M™, (52) 



where 



and 



M 2 J = No^ ^T' -u\ cos(2ttI^) e~ 2 ^ (53) 

7r ti I sinh A; oj c ui c 



2d _ i- I \ / [i H c2 - H l 2 Xi+2 sinh A; 



i /47rr\ nrj 

\ w c / V w c -ff C 2o (I + 2)A; sinh A ;+2 
cos(2^(Z + 2)^) 



with 



A; = 



cos(27r/^f) 
2n 2 lT 



(54) 



All the previous calculations can be done in the same way for a quasi-two dimensional 
system. We give only the final result 



MT = M 2 f. (55) 



The field interval (J51|) in which the results (52), (p3J), (M) are valid is tiny for an ordinary 



type-II superconductor. However for those particular ultraclean materials with very high 
H C 2 and very small £p, where dHvA effect in the superconducting mixed state have been 
observed (see introduction), the value of \f n/u> C 2 is of the order of ten and the theory 
presented here has noticable region of applicability below the upper critical field. 

As the final remarque, we have to remember that /i is the chemical potential less the 
energy band £(, so that the limit (|5l|) is different for different bands and less restrictive for 
the bands with smallest value of fi — Eb ~ £f — £b- 
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5 Conclusion 



The grand canonical ensemble theory of de Haas- van Alphen effect in the normal and super- 
conducting states of 2D and quasi 2D metals valid at fi ^> lo c is developed taking into consid- 
eration finite temperature, level broadening due to impurity scattering, spin paramagnetic 
splitting and additional suppression due to inhomogeneous order parameter distribution in 
the mixed state. The main results are represented by the expressions (p3|), ( |34| ) and (p2|), 
(p33|), ( |5~4| ) for the magnetization amplitude oscillations. 

The straightforward generalization to the multiband case is obtained by the summation 
of these expressions over several energetic bands. 

The limits of observation of dHvA effect near the upper critical field ( |5l| ) in the super- 
conducting mixed state will be of course different for different energetic bands being less 
restrictive for the bands with smallest orbit area ~ {ef — £&)• 
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Appendix A : Calculation of a 

In this Appendix, we perform the calculation of the coefficient a(H, T) in the quasi-classical 
limit. Effects of quantization on the upper critical field have already been studied by Gruen- 
berg and Gunther (1968) for a 3D system. These authors give a method for the calculation 
of the corrections in the linearized gap equation, which can be adapted to our 2D case. 
According to equation ([36]) we have 

+00 

a{H,T) = --TY$a, v (56) 

where 

S av = 23? y e-^G ff (R,c^)G- CT (R,-c^)dR. (57) 

In the 2D case, the Green's function G CT (R, Q v ) in the representation R, after the integration 
over orbit centers (Dworin 1966), is expressed as 



g*(r, a v ) = iW'/s g — (58) 

^ iu) v - £ r + a[i e tl 

where L r (t) is the Laguerre's polynomial of order r, t = Hp 2 /2, and 

Cr = {r + -)ui c - (J,. 
Using the relation between Laguerre's polynomials 

+oc /1\ n+m+1 (m + nV 



" 



wc obtain 



2/ m\n\ 



(r + iy. (l/2) r+l 



fl^O ( lUJv + - VeH) + 6 + Me-n ) 
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Following Gruenberg and Gunther (1968), we do the Gaussian approximation 

1V+' (r + l)\ e-^i^ 
I Hi! M (7rr) V2 

valid for fi/uj c 3> 1. Then we make use of the Poisson's formula which gives us three kinds 
of terms : 

+00 

sa„= E 5 sr> ( 60 ) 

n,m— — oo 

<™ m „ (_ 1) n+m JV gj / / rf * « _ ^ (61) 

JO Jo V 7 ™ + 6; - Me-ff ) +£y + HeH) 

with 

^ x = XLU C - (A. 

In the first term obtained for n = m — 0, we transform from integration over x and y 
to integration over the coordinates of the two-dimensional vectors q and q', so that u> c x — 
q 2 /2m, ui c y = q' 2 /2m. Since q and q' are very close (and of order of ]cf), we have 

{x-y) 2 /Ax^{q-q') 2 /2H 

After using the following approximative relation 

e -(q-q') 2 /2H 1 f-2-K 



^qq'/2H V^Jo 

changing variable q' to Q = q q' and neglecting the term Q 2 /2m in the energy wc 
obtain 

00 Nqcj c r -Q 2 /2H f dq 

Since q is very close to /cp, the term q/m is approximate by the Fermi velocity Vf- The 
integration over q is changed to an integral over energy and angle 9. The first integration 
over energy gives 

S™ = jj^b R [ dQe -Q 2 /2H f 2 " dB 

w » ttU J J 2u) v + 2iLi e H -iv F Qcos9' y ' 

The second integration over the angle is performed by introducing the complex variable 
z = e ld . The corresponding contour is the unity circle. We find 

_, 00 2N ^ f e-Q 2 / 2H dQ 



s m> = r^sjj , ( 6 4) 



S?° = 27rJV CR / (65) 
Jo 



Rewriting this expression as 

r+ °° e~< x dx 

y/n + X 

where ( = 1/sf^c (sf is the Fermi level) and k = (u>„ + i[i e H) 2 , then we calculate 

rp +°° 

coo _ 1 coo 
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Using the following relation valid at low temperatures 



r+oo 2 ff)F\ 



(66) 



we get 



' S 00 = ^- f^lnxe^dx + ^-T 2 ^ [ + °° - 6 - hA (67) 



Here kq = T + i/i e H. In the first term in the right-hand side of equation (|66j), we neglect 
r and fi e H. To avoid the divergence, a cut-off loq is introduced ; this latter vanishes with 
subtraction of the term = ln( 2 Jy° ) where C = In 7 is the Euler constant. Taking notes 
of the value of the integrals (Gradshteyn and Ryzhik 1980) : 



H-00 

C / Inxe'^dx = -C-lnC 
./o 



and 

2V7rCe KoC 1 = / e~ l dt 



Jo (x + k ) 3/2 ^/ko~ " s I V™ Jo J ^f^^o 

we obtain 

i_ T g sS = ^{,„ ( JL ) + ^) (68) 

5 ^0 c2 ° ^ Wc ^ 

with lf C 2o = 2?r a"" ■ The temperature dependent term is neglected further in so far as we 
work at very low-temperatures. 

For the second term obtained for n = m =/= in (|6^), we follow the steps leading to 
equation (|6^). We approximate the additional term in the exponent by 2z7rn(x — y) « 
Imnvp cos 0. Thus we find for n > 0, S~: n = and for n < 



5™ = 2^0^-^'^^ / + °° (69) 



J \/K + X 

Using the value of the integral (Gradshteyn and Ryzhik 1980) 



yjn + x VC \ V™ Jo / V C 

for £«; -C 1, we get 

g r= 2 ^ j R e -^^4^, (70) 

The third and last term in ( |60| ) containing the summation over n 7^ m is evaluated only 
for nonzero values of n and m because S£® and 5£™ are negligibly small. We approximate 
(7rx) -1 / 2 exp[— (x— y) 2 /4x] by (lUc/it^) 1 / 2 in order to have two uncoupled integrations. These 
are nonzero only for n and m negative and yield for a pair (n, m) 

S nra + S mn = , m+m^^jj e -2n\n+m\ (2*±£«£) ,3 r _ TO )A), ( 7 1) 

Reassembling all these different contributions, we have 
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where 



2n 2 T 2 

So = ^ps (73) 
6/j,uj c 



„ +00 +00 

Sl = ^S^^r^^', (74) 
c n =i i/=o 



-Uqq _|_qq _|_qq 

5 2 = M.sr£ £ ^(-l)»+™cos(2 7 r(n-m)^) e - 2 ^ +m »^ ii ). (75) 

n— 1 m>n zy— 

The term Si is the same as in the 3D case (Mineev 2000). We also do the same transformation 
and include it in the upper critical field at low temperatures averaged over the oscillations. 
Like in (Mineev 2000), we simplify the term S2 in (fzHJ) changing summation variables to n 
and m — n = I, and keeping reasonably only the first term in the sum over n 

irT±^ p -2^(i+2)^ ., tt 

S 2 = — V n-; cos(2tt(/ + 2)^1) cos(27r/ii). (76) 

uj c ^ smhA z+2 u c u c 

This is the term which is taken into account in a osc . 

In the quasi 2D case, we can follow exactly the same previous proceeding to calculate 
the coefficient a. The Green's function is expressed in a mixt representation (x,y,k z ) 

G CT (x, y, k z ,u> u ) = AW c e-'/2 g _ (77 ) 

^— ^ — (r + jJCJc + 2t cos(fc z d) + ;U + a\i e tl 

In equation (|57j), in addition to the integration over the two-dimensional vector R there is 
an additional k z integration. At last, we obtain 

- 2 + (I) 2 * - *Vf* - (78) 



where 



5, 



q2d 




C2d 

d ' 

5l ~~d~' 



~ +00 _27r(i+2)^ tt n, 

S f d = — Yi-l f - - cos(27r(/ + 2)^=-) J (2tt/— ) cos(27rf-^). 
w c smhAi+2 w c w c uj c 



Therefore, the oscillating part of a includes an additional Bessel function factor analogous 
to that in quasi 2D normal metal expressions. 



Appendix B : Calculation of (3 

The calculation of the coefficient (3 is performed in the magnetic sublattices' representation 
(Vavilov and Mineev 1997) 

xG°{U>,Z»)K£> (79) 
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where F"™^, = / ^^/ nm (q)/* m , (q)/„- m (q)/*, m , (q) are the matrix elements of the Abrikosov's 
function f. Keeping only the main contribution to the expression ([79]) given by the diagonal 
terms n — n' — to — to', and since F™" w H/(2ir) 2 n (Vavilov and Mineev 1997), we have 

m T ) = ^E r EE ' yK " a>) ? B, " ,)ia (so) 

^ ' cr±l i/ n=0 

The use of the Poisson's summation formula to transform the summation over n leads to 
a smooth non-oscillating field dependent part and to a fast oscillating field dependent part 
of/3. The smooth part is obtained by the substitution of the summation by the integration 
with 



and n ~ fj,/u> c '■ 



~ P{H,T) 47TA* S T ?/oo [^-e + a Me i7)HcD,-^-a Me i7)] 2 - 



<T±1 U 

The integration over energy £ gives 



This expression is then performed at T = K to allow the change of the summation to the 
integration. We find 



ffi°)^o^ ^:^r (83) 



^ r 2 - ( Me g) 

' 167T/1 (r 2 + (/Z e iTj 

The oscillating part of (3(H, T) is given by 



" CT±1 !/ i=l 

+°° e ~ 27 ™'t^: 



(iOv - £ + (JH e H) 2 {-iuj u - £ - G[i e H) 2 ^ 
The integration over £ presents no difficulty 

,, +oo +oo 

/W#, T) = ^± V(-l) z cos^Z^T V e -™(*^) 

x \ T-~ : TrvT + — 7T : ^r-r ^ • 85 

At T = 0, the summation over is replaced by the integration, which yields at last 

0"(H,T = 0) = ^_g(-l)'cos(2 7 r/A) /(27r/ il) e --^ (86) 



where 

I{x) • ' 



{y + 1)3 (l + y ) 2 
For a quasi 2D system, the same calculations yields 



-xy 



dy. 



^ 2rf (if,0) = ^ 2d( f' 0) (87) 
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and like in the oscillating part of a, an additional Besscl's function factor appears the 
oscillating part of (3 

= T2^ B-l)'-P<)W^»(^)e-*<. (88) 
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